Spectral universality of strong shocks accelerating charged particles 
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As a rule, the shock compression controls the spectrum of 
diffusively accelerated particles. We argue that this is not so 
if the backreaction of these particles on the shock structure 
is significant. We present a self-similar solution in which the 
accelerated particles change the flow structure near the shock 
so strongly that the total shock compression may become ar- 
bitrarily large. Despite this, the energy spectrum behind the 
shock is close to E~ 3 ^ 2 independently of anything at all. 
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The first order Fermi acceleration at shocks was al- 
ways regarded as one of those universal phenomena that 
result in a scale-invariant behavior of physical quantities. 
In particular, this acceleration mechanism has predicted 
a power-law momentum distribution ^,|^] similar to those 
observed in diverse astrophysical objects. There exists, 
however, another important characteristics, the energy 
content of accelerated particles. Generally speaking, this 
quantity makes the expectations of the scale invariancc 
not so well grounded. Whenever the shock accelerates 
efficiently in the sense of its energy conversion into ener- 
getic particles (and this is precisely what one anticipates 
in the context of the problem of cosmic ray origin), a 
more complicated flow structure must form in place of 
this shock. Besides the shock itself, usually called sub- 
shock, it contains a long precursor in which the inflow- 
ing gas is gradually decelerated by the pressure gradient 
of energetic particles, in other words of the cosmic rays 
(CRs) produced in the entire shock transition. Fig. 1. 
One interesting aspect of this gas-CR coupling is an in- 
crease of the total compression ratio which should harden 
the spectrum and make the CR acceleration potentially 
much more efficient. At the same time, the most attrac- 
tive feature of this acceleration mechanism, the scale- free 
form of the spectrum is likely to be lost since the shock 
develops now its internal structure and particles with dif- 
ferent momenta "see" different shock compression due to 
the strong momentum dependence of the CR diffusivity. 
It is also well known (see e.g., Ref. 0) that this broad- 
ening of the shock transition softens the spectrum. The 
combined effect of these two oppositely acting factors has 
been resisting any regular analytic calculations for a long 
time (see, however, || for a perturbative approach). 

Recently Q , we have reduced the problem of the gas- 
CR coupling to one nonlinear integral equation. Its solu- 
tion does provide selfconsistcntly both the particle spec- 
trum and the structure of the hydrodynamic flow. At 
least in a certain phase space region this solution turns 



out to preserve the scale-invariance of the test particle 
results. But the most surprising aspect of this solu- 
tion is that in contrast to the well known linear solu- 
tion for the particle momentum distribution / oc p~ q 
where q = 3r/(r — 1) and r stands for the shock compres- 
sion ratio, in the nonlinear regime q = 3^ independently 
of anything at all. This means that when accelerated 
particles modify the shock strongly their spectrum and 
the flow profile mutually adjust in such a way that all 
such shocks become universal accelerators producing es- 
sentially the same power-law spectra. They differ only in 
the cut-off momentum that is always limited by the size 
and/or the age of the shock. The above result is, however, 
limited formally to relativistic particles with the Bohm 
diffusivity n oc p and relates to the solution (out of the 
three possible) with the highest conversion efficiency of 
the flow energy to the CR energy. 




FIG. 1. The flow structure in a strongly modified CR shock. 

In this letter, using a simpler mathematical tech- 
nique, we demonstrate that the spectrum universality 
persists also in the nonrelativistic region and for even 
more general form of diffusion coefficient, e.g., n oc 
p 2 1 yl + p 2 /m 2 c 2 . The power-law index, being q = 3^ in 
the region p 3> mc, simply changes for q = 4 in the region 
p <C mc before it conforms to the local subshock compres- 
sion r s , i.e. approaches q = 3r s / (r s — 1) at p ~ po <C mc, 
where po is the injection momentum |^|,|| . Remarkably, 
the q — 3^ index coincides with the plain test particle 
result for a strong shock in a purely thermal relativis- 
tic gas, whereas q = 4 coincides with that for a non- 
relativistic gas. All these in spite of the facts that the 
flow profile is strongly modified and the total compres- 
sion ratio may be much higher than 4 (q = 4) and even 
7 (q = 3^) occurring in nonrelativistic and relativistic 
gases, respectively. Moreover, the subshock compression 
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ratio r s may be significantly lower than 4 and cannot 
be responsible for q = 4 spectrum in the nonrelativistic 
momentum range. 

The standard formulation of the problem includes the 
diffusion-convection equation constrained by conserva- 
tion of the fluxes of mass and momentum (see e.g., 
Rcf. 0). We assume that a strong CR- modified shock 
propagates in the positive x-direction in a cold medium 
with the mass density p\. In the shock frame of refer- 
ence the steady mass flow profile is defined as follows 
U(x) = —u(x), x > and U(x) = —U2, x < 0, where 112 
is the (constant) downstream mass velocity, u(0+) = u , 
and u(x) — > u± = const, as x — > 00, Fig. 1. The equa- 
tions read 
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dg\ 1 du dg 



pu = p x u\, 
P c + pu 2 = p x u\, 



x>0 



(1) 

(2) 
(3) 



Here the number density of CRs is normalized to 
47r gdp/ p. The upper limit pi stands for a boundary 
in the momentum space (cut-off) beyond which parti- 
cles are assumed to leave the system instantaneously 
(g = 0, p > pi). In the downstream medium, x < 0, 
the only bounded solution is g — G(p) = g(p, x — 0). In 
Eqs. (0,||) p(x) is the mass density, p\ — p(oo), P c is the 
CR pressure 
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(4) 



The particle momentum p is normalized to mc. As in- 
dicated, Eq. (|^) is written in the region x > where 
we have neglected the contribution of the adiabatically 
compressed cold gas (i.e., particles with < p < p ) con- 
fining our consideration to sufficiently strong shocks with 
M 2 = p\u\/^P g i ^> (ui/ito) 7 , where 7 is the specific 
heat ratio of the plasma and P g i is the gas kinetic pres- 
sure far upstream, i.e., at x = 00 (sec Q] for a detailed 
discussion of this approximation). The subshock strength 
can be obtained from the familiar Rankinc-Hugoniot con- 
dition for the gas 
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7 + 1 



u 2 7- l + 2(ui/w ) 7+1 M- 2 



(5) 



Introducing the flow potential 0, such as u — d<j)/dx 
we seek the solution of Eq. (|l|) in the form 



9 = 5o (p) exp 



1 + /? 
n(p) 



x>0 



where 



(3{p) = -{l/3)d\ng /d\np 



(6) 



(7) 



This is the key substitution in our analysis. Considering 
u as u(4>), substituting Eq.(|^) in Eq. ([!]) and separating 
the variables we obtain 



^ = Xu 



df] 

P-T- = (1 

dp 



dhiK 
dlnp 



(8) 



(9) 



where A is a constant. Eq. (g) may be readily integrated 
and yields for the flow potential 



<P(x) = A/(1 ~ A) [(1 - X)u x + M 1/(1 ~ X) 



(10) 



where (j> = (f)(0) is another constant (see below). 

The solution for the particle spectrum go(p) to be ob- 
tained from the system (0,^|) depends on two further inte- 
gration constants. These are the magnitude and the slope 
of go at the injection momentum p — po, i.e., go(po) and 
fto = (3(po), respectively, provided by an "injection" solu- 
tion 0. Throughout this paper we shall consider goipo) 
as given. The constant Po may also be determined from 
merely the subshock strength (see also ref. B and be- 
low). It is straightforward to verify that the following 
expression is the first integral of the system MM) 



9o(p) 



1 + /3 



= 3o (po) 



K>0 



1+00 



(11) 



where kq = n(po)- The solution for go is given by 



9o(p) = 9o bo) 
Po + 
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Xk 



Po 3/X I K ip')P 

Po 



1 -A 



,/3/A-l 



dp' 
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For p not too close to po, more precisely for 
(k/ko)(p/po) 3 ^ X ^ I, the particle spectrum g acquires a 
form determined by merely that of n(p) and "forgets" its 
behavior at p ~ po- Note that the latter (and thus the 
parameter Po) is prescribed by the subshock strength, 
according to the relation d\nG/d\np = 3/(r B — 1) at 
p = po , which can be recast with the help of (nAo as 



Po - 4>oPo 



l + Po 

K 



1 
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(13) 



If p is large in the above sense and k(j>) is a power-law, 
then the spectrum go(p) must also be, go{p) oc n~ x (p) 
and 



P~ (A/3)dln K /(ilnp 



(14) 



(the last formula may be also obtained from Eq.(||) as- 
suming P ~ const). As we shall see, the parameter A 
depends on the scaling of k as well, and the most surpris- 
ing consequence of this dependence is that the resulting 
slope of go{p) is, in fact, independent of k(j>). 

What we obtained so far is a one parameter family of 
exact solutions to equation (Q) that require a rather spe- 
cial class of the flow profiles u(ifi). It is by no means guar- 
anteed that any of these solutions satisfy the Bernoulli's 
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integral (||) (the continuity condition (||) can be easily 
satisfied since it does not depend on g(p > po)). Indeed, 
the only unspecified quantity in solution (j|,|l^,|l^) that 
can potentially be used to satisfy the functional relation 
(|J), is the constant A which is, generally speaking, not 
enough. Fortunately, the dependence of P c {4>) in Eq.(||) 
can be made consistent with the dependence u(<f>) oc <^> , 
following from Eq. (|8|) , in a fairly large part of shock tran- 
sition, provided that the parameter A is chosen properly. 
To demonstrate this we substitute Eq.(p|) into Eqs.(||,|). 
Using Eq.(0), the Bernoulli's integral (0) rewrites 



»(0) + c 



P 2 {s) 



VTTpHs) 



Ul 



(15) 



where we have introduced a new variable s in place of p 

s=^ (16) 

K 

and the limits sq,i = s(po t ±). We have also used the first 



integral (|1 
here as 



, go oc s . The injection rate £ is defined 
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to 
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(17) 



and the function p(s) in Eq.jl^) should be determined 
from Eq.©. We now make a "standard" assumption: 
n{p) = Kp 2 (l +p 2 ) -1 / 2 , i.e., the mean free path of a 
particle is proportional to its Larmor radius (here K is a 
reference diffusivity) . Then, Eq.([l5|) rewrites 



u((j)) = ui 



K 



S x ~ 2 e-t> s ds 



(18) 



According to Eqs.( p^ , ^6| ), (3(s) is a very simple function, 
taking in the most part of its domain nearly constant 
and relatively close values, f3 ~ A/3 for Ks <C 1 and 
(3 ~ 2A/3 for Ks 3> 1. It changes monotonically between 
these limiting values. Therefore, the last integral allows 
a straightforward asymptotic analysis. Our goal here is 
to examine whether Eq.(||), according to which du/dtfi oc 
A_1 is compatible with Eq.(|l8|). We assume that < 
A < 1, which will be confirmed later. Differentiating 
Eq. ( |l8| ) with respect to c6 and considering first the region 
1/so <C <t> <C 1/si we may obviously replace the lower 
limit by zero and the upper one by infinity. From Eq. ( |l6| ) 
we then obtain 



du 



Kcj) x 



[l + l/l3(T/<P)]T A - l e- T dT 



§p [i + mm] 



(19) 



where T is the gamma function and /3(r/c6) is replaced 
in the last integral by its mean value at f/<p with f ~ 1. 
As we have seen already the function /3(f /</>) in the last 



expression varies slowly and it is close to A/ 3 for <f> > K 
and to 2A/3 for <f> < K. Therefore, the 4> dependence 
of du/d(f> is determined by the factor <f>~ x and is indeed 
consistent with Eq.(||), i.e., with u oc A for A = 1/2. 
Eq.(|l9|) becomes invalid for (j) l/ s i> since the lower 
limit in Eq.(|l8|) cannot be replaces by zero in this case. 
The function du/d<p cuts off as (see Eq.(nS)) 



du 



-si<t> 



(20) 



The last asymptotics obviously matches the formula ( |l9| ) 
at (j) ~ 1/si but it is no longer consistent with Eq.(||). At 
the same time, this is a periphery of the shock transition 
where u((f>) exponentially approaches its limit u% at x ~ 
I = n(pi)/ui, Fig. 1. This particular form of the cut- 
off of du/d(j) is entirely due to our assumption about an 
abrupt cut-off of the particle spectrum at p — p\ (s = si 
in Eq.(|l8|)). Therefore, there is no major physical reason 
to refine our solution in this region unless a more realistic 
model of particle losses in the region p ~ p\ is adopted. 
According to Eq. (|l0|) the flow profile in the internal part 
of the precursor is simply linear (A = 1/2) 



U(X) =Uq + ——X 

20o 



(21) 



An asymptotic approach resulting in a flow profile u(x) 
that is uniformly valid (for all x > 0) and coincides with 
cq.(|2l]) for < x < I has been developed in Ref. [Q. 
It produces a full set of Rankine-Hugoniot relations, de- 
termining the total (ui/112) and the subshock (1*0/^2) 
compression ratios along with the constant <f>$ in Eq. ( |l0| ) 
as functions of the Mach number, the injection rate and 
the upper cut-off momentum p\. These quantities have 
been calculated in Ref. Q| with a help of nonlinear inte- 
gral equation derived from the system (|l|j^). This equa- 
tion has three different solutions in a certain parameter 
range (see also 0) with dramatically different acceler- 
ation efficiencies 0] . However, its derivation utilizes a 
relativistic form k oc p for all p and a similar form of 
the pressure integral. These assumptions oversimplify 
the particle spectrum at p < 1 and we have dropped 
them here. According to the form of the particle spec- 
trum with the index q = 3(/3+ 1) = 3.5 at high momenta 
p > 1 and the linear profile of the flow velocity ( pl| ) the 
solution obtained in the present paper can be identified 
with that of the Ref. Q having the highest efficiency of 
CR production. Note that this solution becomes unique 
for sufficiently large ( and pi . In addition, it provides pa- 
rameters that have not been determined in the present 
treatment. First of all the relation between the flow de- 
celeration in the precursor and in the subshock was found 
to be (ui — Uq)/[uq — 1*2) oc Qp\ which, in combination 
with Eq. (||) determines the unknown velocities Uq and u-i 
given the far upstream velocity u\ . Also the constant (f> 
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may be related to the precursor length I — k(p\)/ui by 
ul/2cj) ~ ui/l. 

As we already emphasized the scale-invariance does 
not exist everywhere in the phase space as in any other 
solution of this kind applied to a real physical situa- 
tion. However, the scale invariant behavior in the re- 
gion Mo < u < ui, pa < p < pi may span decades in 
both momentum and coordinate space in those astro- 
physical situations where the dynamical range pi/pa is 
large. E.g., for typical supernova shock conditions one 
may expect px/po ~ 10 8 - It is also interesting to note 
that the downstream particle spectrum given by Eq. ( p"l] ) 
with the approximate formula (|l4|) for j3 and being ex- 
pressed in terms of kinetic energy E rather than momen- 
tum exhibits a particularly uniform behaviour through- 
out the entire energy range, relativistic and nonrelativis- 
tic. In a standard normalization F{E)dE this spectrum 
has the form 



(E + 1)(7E 2 + UE- 



(£ + 2)3 



where E is measured in mc 2 , and the spectral index is 
close to 1.5 for all but injection energies (if the subshock 
is substantially reduced, see text below Eq.(|l2|)). We em- 
phasize that its value does not depend on any parameters 
involved in these calculations. This results from the mu- 
tual adjustment of the flow profile u(x) and the particle 
distribution in energy and in this sense the form of the 
spectrum is completely universal. 

One may ask then what does the power-law index de- 
pend on? It should be primarily the momentum depen- 
dence of the CR diffusivity n(p) that we have specified 
in our treatment above. To examine this idea we rescale 
k as follows k' = K a , where a > 1/2 and the rescal- 
ing of all primed variables below is induced by the above 
transformation of k. According to Eq.(|l4|) the spectral 
slope (3 is now to be replaced by (3' — (a\' /3)dlnK,/d\np. 
Recalculating du/d(f> in Eq.(|l9|) with these rescaled spec- 
trum and CR diffusivity k' we obtain du/d<j) oc c^ 1 /" -1- * . 
Since the formula u cx <j) X (Eq.(||)) holds, we deduce that 
the rescaled A' = l/2a. Therefore the spectral slope j3 
remains actually unchanged, f3' — [3. 

We conclude that the spectral universality is insensi- 
tive to the scaling of the CR diffusivity «(p), in other 
words, to the spectrum of the underlying MHD turbu- 
lence. On the other hand, the velocity profile does de- 
pend on a. From the above analysis we obtain u oc 
a .l/(2a!-l) which a i so explains the condition a > 1/2 ||. 
There is a small deviation from this scaling at the dis- 
tances x corresponding to the diffusion length of particles 
with p ~ 1 since (3 depends on </> in Eq.© at K. 

Summarizing the last results, when n{p) rescales, so 
does the flow profile u{x) but (3 remains invariant. It is 
not difficult to understand why this is so. As usual in the 
Fermi process the spectral slope of course depends on the 



flow compression. But, since the flow is modified, a parti- 
cle with momentum p, samples not the total compression 
but only a compression accessible to it. This compression 
is determined by the relation <p(x) oc n a (Eq. @). As 
we have shown, u{4>) oc Therefore, the flow com- 

pression, as seen by this particle diffusively bound to the 
shock front, scales as u/u2 oc ^l 2a oc yjn{p). As this is 
independent of a the index (3 must also be. 

This work was done within the Sonderforschungsbere- 
ich 328 of the Deutsche Forschungsgemeinschaft (DFG). 
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